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Abstract 

A key element in the generation of optical torque in optical traps, which oc- 
curs when electromagnetic angular momentum is transferred from the trap- 
ping beam to the trapped particle by scattering, is the symmetries of the 
scattering particle and the trapping beam. We discuss the effect of such 
symmetries on the generation and measurement of optical torque in optical 
tweezers, and some consequent general principles for the design of optically- 
driven micromachines. 
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1. Introduction 

Optical t rapping, exernplified by the single-beam gradient trap, or opti- 



cal tweezers ( lAshkin et all Il986l ). has developed from an interesting novelty 



into a widely-use d and powerful tool with n i any a pplications in biophysics 



and other areas (lAshki: 



1 p owern 
M I2OOOI : 



Moffitt et al.l . l2008l ). Uses include the mea- 



surement of forces on the order of picon ewtons, such as the adhesion forces 



on a cellular level (Knoner et al 



20061 ) and forces in single-molecule stud- 



ies ( IWilliams and Rouzinal . |2002| ). as well as manipulation of microscopic 
objects, including live bacteria and eukaryote cells. 

An optical tweezers trap is conceptually simple: a high numerical aperture 
lens, typically a microscope objective lens, is used to focus a laser beam to 



Email address: timo@physics.uq.edu.au (Timo A. Nieminen) 



Preprint submitted to JQSRT 



December 11, 2008 



a diffraction-limited spot. The sample, consisting of the microscopic objects 
one intends to trap or manipulate (often polystyrene or silica microspheres) 
suspended in a fluid, which would usually be water or a biological buffer 
solution, is placed on a stage, and the objects, as long as their refractive index 
is higher t han that of the susp ending fluid, and they are not too reflective or 



absorbing (jStilgoe et al.l . l2008l ). can be trapped in the bright focal spot. 



The optical forces responsible for trapping are fundamentally the result 
of the trans fer of momentum from the beam through scattering of the trap- 



ping beam (lAshkin et al.l . Il986l : lAshkiru . 1 19921 : iNieminen et al.l . l2007al ) ; this 



possible due to the transport of momentum by electromagnetic fields, with 
the optical force being equal to the rate of transfer of momentum from the 
beam to the particle. 

Since angular momentum can also be carried by electromagnetic fields, in- 
cluding light, the possibility of exerting optical torque on particles, trapped 
or otherwise, is evident. This introduces the possibility of six-degrees-of- 
freedom, with the orientation of a trapped particle being controlled, in addi- 
tion to its position — this would be a major advance in the manipulation possi- 
ble with optical tweezers. Rotation can also be used as a tool to actively probe 



Bishop et al.l . 12004 iLa Porta and Wand . |200J). In addition, i t could allow 



m 

the physical properties of the surrounding me dium (INieminen et al.l . 12001 



the development of practical optically-driven microma chines ( iFriese et al. 



200ll : iGalajda and Ormosl . l200ll : INieminen et al.l . l2006a ) . A significant amount 



of progress h a s been achieved in all of the se areas (e.g. 



Parkin et al.l . l2007l : iKelemen et al.l . |2006| ). using a diverse range of methods 



Bavoudh et all . 12003 



to generate the optical torque. 

A key element determining whether or not the transfer of electromagnetic 
angular momentum to a particle by scattering is possible is the symmetries 
of the scattering particle and the trapping beam (iKonz and Benfordl . 12003 ; 



Nieminen. 2004 



Benford and Kong . 12004 : INieminen et al.l . 12004 ). Therefore, 



consideration of the effects of symmetry can provide insight into rotational 
optical micromanipulation and design of optically-driven micromachines and 
a unified overview of the generation of optical torque. 



2. Optical angular momentum 



Optical forces have a long and interesting history (jWorralll . Il982l ). start- 
ing with Kepler's conject ure tha t radia tion pressure was responsible to the 
direction of comet tails (iKeplerl . Il619l ). However, due to their smallness. 
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practical applications in the laboratory did not appear feasible until Ashkin 
noted that only ver y small forces are needed to manipulate very small ob- 
i ects (lAshkinl. |2000[) . which, notably, led to the invention of optical tweezers 



(lAshkin et all Il986l ). The fundamental principle by which optical tweezers 



operate is the ex change of momentuni between the trap beam and a particle 



within the trap (lAshkin et al.l . Il986l : iNieminen et al.l. l2007br) 



It was suggested as early as 1909 by iPoyntind (119091 ) that light can carry 
angular momentum as linear momentum. While Poynting was skeptical as 
to whethe r such tor ques c o uld b e measured, this was achieved in 1936 by 



HolbournI fll936t ) andlBeth] fll936[ ). 



Most of the earhest optical tweezers-based demonstrations of optical torque 
made use of absorption, combined with la s er beams carryi n g ang ular momen- 



tum (IHe et al.l . Il995l : iFriese et al.l . Il996l : ISimpson et al.l . 119971 ). which has 



very limited practical use due to the resulting heating, which, for example, 
could lead to rapid death of live specimens, and in any case, makes three- 
dimensional trapping more difficult or even impossible. 



More pronaising was the use of birefri i igent materials (iFriese et al.l . 11998 



Bishop et al.l. 12004 iLa Porta and Wang. 



(iHigura shi et al 
[2003; Neale et all 



19941. ll997l : lBonin et al. 



20041) , or specially-shaped objects 



20021 : ICheng et all . l2002l : iBishop et al. 



20051 ) which c an themselves be ma de of birefringent ma- 



terials, or be shape-birefringent (IBorn and Woli 119971 ). 

Notably, these types of methods, using birefringence (due to either the 
material or shape) which affects the polarization of the trapping beam and 
therefore depends on the spin angular momentum of light, and the use of 
shapes which depend on, and affect, the spatial structure of the trapping 
beam, and therefore make use of the orbital angular momentum of light, are 
often viewed as distinct physical mechanisms. 

It is true that there are important differences between spin and orbital 
angular momenta, but few of them will be of importance here. One dif- 
ference that will turn out to be important is that spin angular momen- 
tum can be readily determiri ed by measurement of the Stokes parameters 



fICrichton and Marstonl . l2000f ). However, it is still worth reviewing the ba- 
sic elements of electromagnetic angular momentum, especially since there 
is ongoing debate within the literature conce r ning the angular momentuni 



of cir cularly polarized waves (iHumbletl . Il943l : IStewartl . l2005l : iPfeifer et al. 



20061 ). 



The controversy arises since there are two paths we can take to deter- 
mining the angular momentum density or flux. Firstly, we can assume that 
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the angular momentum density is the moment of the momentum density, 
r X (E X H)/c, where r is the position vector relative to the origin of our 
chosen coordinate system, and E and H are the electric and magnetic fields 
respectively, which leads to the curious result that a circularly polarized 
plane wave carries zero angular momentum in the direction of propagation, 
in contradiction to the usual quantum mechanical result of ±h angular mo- 
mentum per photon. This choice of angular momentum density implicitly 
states that there is no spin angular momentum — the spin angular momentum 
is the part of the total angular momentum which is independ of the choice of 
origin about which moments are taken; taking r = in the above expression 
for the angular momentum density yields a spin of zero. 

Secondly, we can begin with the L agrangian for the electromagnetic ra- 
diation field, and, via Noether's (119181 ) theorem, obtain the canonical stress 
tensor and an angular momentum tensor. The angular momentum tensor 
can be divid ed into spin and orbital components (jJauch and Rohrlichl . Il976 ; 



Soperl . Il976l ) , based on the independence of the spin from the choice of ori- 



gin. For a circularly polarized wave, the spin would be I/w, where / is the 
irradiance and u the angular frequency, in disagreement with the first result, 
and in agreement with the quantum mechanical result. This separation of 
the angular momentum density into spin and orbital parts is, in general, not 
gauge-invariant, and it is common to transform the canonical stress tensor 
into a summetric stress tensor, with a gauge-invariant r esult for the angular 
mome ntum density, yielding the integral of rx (ExH)/c. iJauch and Rohrlich 
( 119761 ) carefully point out that this transformation requires the dropping of 
surface terms at infinity. The reverse of this procedure, obtaining the spin 
and orbital terms starting from r x (E x H ) /c, in volving the same surface 
terms, had already been shown by iHumbletl (119431 ). For and physically re- 
alisable beam, these surface terms vanish, and despite the controversy, both 
expressions for the angular momentum density yield the same total angular 
momentum and torque exerted on an object (jZambrini and Barnettl . 12005 : 



Nieminen et al.l . l2007d ). 



Since we are interested in torque within optical tweezers, where we have 
a monochromatic wave, the gauge-invariance of the separation of the an- 
gular momentum into spin and orbital parts is also a non-issue, s ir ice in 
this cas e , the separation is gauge-invariant ( Crichton and Marston . 2000l : 



Barnettl, [2002) 



Consequently, the spin angular momentum flux of a paraxial electromag- 
netic beam, where the wavefronts are almost plane (by definition), varies 



4 



from —h per photon in the direction of propagation for right circular polar- 
ization, through zero for linear polarization, to +h per photon for left circular 
polarization. In classical terms, these extremes are ±P/uj, where P is the 
beam power. Linear or elliptical polarization can be represented by the su- 
perposition of two circularly polarized components of opposite handedness; 
the total spin angular momentum is equal to the sum of the momentum of 
the two components, which, if they are of equal power, will be zero. 

Since, in the paraxial limit, the spatial structure of the beam is inde- 
pendent of the polarization (since it is described by a solution of the scalar 
paraxial wave equation) , the spin and or b ital a ngular momenta are inde- 



pendent of each other (j Allen et al.l . 1 19921 . 119991 ). Thus, an excellent for 



malism for the description of the angular momentum of paraxial beams is 
provided by the combi nation of Laguerre-Gauss beam mod e s and circularly 



polar ized components (iHeckenberg et al.l . Il999l : I Allen et al.l . Il999l : iNienhuisl . 
20081 ) . Each Laguerre-Gauss mode is described by a radial mode index p 



and an azimuthal mode index i, and is denoted LGp£. Each mode has an az- 
imuthal phase dependence of exp(i£0), where is the azimuthal angle. This 
azimuthal dependence is typical of that in solutions to PDEs in rotational 
coordinate systems obtained by separation of variables. Each mode also has 
a well-defined orbital angular momentum flux of iH per photon about the 
beam axis, or iP/u, which is also typical of such solutions. In this case, the 
modes are orthogonal w.r.t. orbital angular momentum, and the total orbital 
angular momentum about the beam axis, L^, can be found by adding the 
individual angular momenta of each mode: 



oo CO 



p=0 ^=-oo 

where Qpi are the mode amplitudes, and units are chosen such that the total 
power is equal to 

CO oo 

|2 



^ = E E M"- (2) 

p=0 l=—oo 

The total angular momentum flux about the beam axis, J^, is then the sum 
of the orbital (L^) and spin {Sz) components, 

J, = + 5„ (3) 

where 

Sz = a,P/uj, (4) 
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where is the degree of circular polarization, which varies from —1 for right 
circular polarization to +1 for left circular. 

However, for a tightly focussed beam such as employed in optical tweezers, 
the paraxial approximation cannot be used — it is necessary to use modes 
which are solutions of the vector Helmholtz equation. It is convenient to 
begin with a set of solutions to the scalar Helmholtz equation, tpni from 
which a set of solutions to the vector Helmholtz equation can be obtained: 

(5) 
(6) 

(7) 



Ln — 




M„ = 


-a X L„, 


N„ = 


1 

-V X M, 

k 



where a is a unit vector or a constant vector (jBrockl . l200ll ). The L„ wavefunc- 
tions are curl-free, while the M„ and N„ wavefunctions are divergence- free. 
Therefore, for electromagnetic fields, we only require M„ and N„. We can 
also note that M„ = {l/k)V x N„. 

In spherical coordinates, we can readily write a general solution to the 
scalar Helmholtz equation: 



(8) 



n=0 m= 



where a„m are the mode amplitudes, j„ are spherical Bessel functions, and 
are normalised spherical harmonics. The resulting vector solutions 
are the regular vector multipole fields, or vector spherical wavefunctions 
(VSWFs): 



i{kr) 
.{kr) 



NnJn{kr)C, 
krN„ "'"^ 



(9) 



+ Nn{j„-i{kr) 



njnjkr) 
kr 



Br 



where = [n{n + 1)] ^^'^ are normalization constants, and B„m(^, 0) = 
rVi;"^(^,0), CnmiO,(l>) = V X ( rYr(e,<i))), and P „, p (g,( / . ) = rF^f^^^) are 



the vector spherical harmonics ( Mishchenko et al. 



20001: Watermar 



Mishchenkol . Il99ll : Ijacksonl . \l99i 

Since the spherical harmonics Ynm have an azimuthal dependence of 
exp(im</)), we once again obtain angular momentum about the beam axis 
in terms of the azimuthal mode index (m); in this case we have the total 
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angular momentum about the beam axis equal to mh per photon, or mP/uj, 
rather than just the orbital angular momentum. The radial mode index n 
also relates to the angular momentum — the VSWFs are eigenf unctions of 
the angular momentum operator J^, with eigenvalues [n{n + 1)]^''^, as well 
as the angular momentum operator J^, with eigenvalues m. As a conse- 
quence of this, it is clear that we must have n > \m\. The radial behaviour 
of each VSWF mode depends on spherical Bessel functions jn{kr), which 
are peaked in the vicinity of kr ^ n. Therefore, the minimum width of 



a focu ssed beam depends on its angular momentu m flux (jNi eminen et al.J, 
2008al ): this has been noted both exp erimentallv (jCurtis and Grie r. 2003j) 
and on elementary theoretical grounds ICourtial and PadgettI (120001 ) . An in- 
teresting result of this, since the minimum width of a beam also depends on 
its wavelength, which also affects the angular momentum flux due to the de- 
pendence of angular momentum on the frequency (or in quantum terms, the 
angular momentum per photon is constant, and the photon flux increases as 
the frequency is decreased) , this minimum width is independent of the wave- 
length and whether the angular momentum is spin, orbital, or a combination 
of the two. 

Since spherical harmonics are a common feature of solutions to PDEs 
obtained by separation of variables in spherical coordinate systems, they, 
and their relationship to angular momentum, are discussed in m any s o urces 
most commonly in the context of quantum mechanics 
Varshalovich et aP . 1988 ). 



e.g., iRosd . Il957 



3. Symmetry, scattering, and angular momentum 

Since we are already considering the trapping be am as a superposition o 
VSWFs, it is natural to use the T- matr i x formalisrn (iNieminen et al.', '2003aj: 



Kahnertl . I2OO3I : iMishchenko et al.l . [200ol : IWatermanl .[l971: Mishchenko . 1991 
Jacksonl . Il999l ). For this, we need to be able to describe the scattered wave, 
for which we use the outgoing VSWFs instead of the regular VSWFs; to 
obtain these, we simply replace the spherical Bessel functions jn{kr) in the 
regular VSWFs with spherical Hankel functions of the first kind hn\kr) 
(assuming a time-dependence of exp(— ia;t) — if we have chosen a time de- 
pendence of exp(ia;t), we use spherical Hankel functions of the second kind 
h^n\kr)). 

Since "T-matrix method" is often used synonymously with "extended 
boundary condition method" (EBCM) which, properly speaking, is a method 
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of calculating the T-matrix (jNieminen et al.l . l2003al : iKahnertl . |2003| ). it is 



worth briefly describing the T-matrix description of scattering. Essentially, 
we expand both the incident and scattered waves in terms of a sufficiently 
complete basis set of functions {ipn^'^^ for the incident wave, and ijjn^^^^ for 
the scatterered wave, where n is a mode index labeUing the discrete modes), 
each of which is a solution of the Helmholtz equation: 



CO 

f/inc = J]a„V'f^), (10) 

n 



t/scat = 5^P.^f"\ (11) 
k 

where the expansion coefficients for the incident wave, and pk are the 

expansion coefficients for the scattered wave. If the response of the scatterer 
is linear, the relationship between the incident and scattered waves can be 
written as a linear system 

CO 

Pk = ^Tknan (12) 

n 

or, in more concise notation, as the matrix equation 

P = TA (13) 

where Tkn are the elements of the T-matrix. 

Note that apart from the assumptions that the scattering is elastic and 
monochromatic, there are no other restrictions other than that a suitable 
basis set exists. This formalism can be used for a wide range geometries, 
including compact particles (ideally treated in spherical coordinates), gener- 
alised cylindrical particles and surfaces, and different types of waves, includ- 
ing scalar and vector waves. 

Since we are interested in the transfer of angular momentum to the scat- 
terer, it is convenient to make use a basis of incoming and outgoing VSWFs, 
rather than regular and outgoing. The incoming VSWFs are simply obtained 

by replacing the spherical Bessel functions jn{kr) in the regular VSWFs with 

(1') 

spherical Hankel functions of the second kind hn {kr). 
Thus, while we began with the incident field written as 

oo n 

Einc(r) = E E ^nlMnmikr) + h^l^^kl), (14) 
n=l m=—n 
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we can replace this with the purely incoming part of the incident field, 

oo n 

Einc(r) = E E ^^l^'nlikr) + b^lN^l{kr). (15) 

n=l m=—n 

In both cases, the scattered field is written 

oo n 

Escat(r) = E E Pi'i.MW (^r) + g« N« (fcr), (16) 

n=l m=—n 

but it is important to note that the "scattered" field is not the same in the 
two cases — in the first case, it is the scattered field in the usual usage of the 
term, while in the second, it is the total outgoing portion of the field, which 
includes the outgoing part of the incident field. Since j„ 
aim = 2aim and bnm = '2bnm, and the purely outgoing mode coefficients for 
the second case (aim""*, ^'■^'''°"*) are related to the scattered field amplitudes 
in the first case {a'nm''^^ ,b^^^'''''''^) by ai^2^°"* = + a'nm'''^^ and = 

&nm/2+6nm**'^^*- The T- matrix itself will also differ in the two cases — using the 
regular expansion, the T-matrix in the absence of a scatterer is a zero matrix, 
while using the incoming field expansion, the no-scatterer T-matrix is the 
identity matrix. The relationship between them is t(™/°"*) = 2T*^"^'^/'^'^'^*^ + 1; 
the T-matrix for the second case is typically written as S and called the 
S-matrix. 

The angular momentum fiux into or out of the system can be found by 
integrating the moment of the momentum fiux density 

j = lRe(r X (E X H*)) (17) 

over a sphere surrounding the scatterer. Taking care to find the moment 

first, and then taking the far-field limit, this integral can be performed ana- 
lytically if the fields are written in terms of VSWFs, for normalised angular 
momentum about the 2;-axis of 

oo n 

= E E + ~ \Pnm\'^ - \qnm\'^)/P (18) 

n=l m=—n 

where 

oo n 

^=E E Kml"" + Kml'' (19) 
n=l m=—n 
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is the incident power. This is the optical torque exerted on the scatterer, in 
units of h per photon. T he spin and orbital con t ributions to the torque ca n 



be calculated if desired (ICrichton and Marstod . l2000l : iBishop et al 
the spin torque is 



20031); 



^ oo n 

4ee 



m 



P ^ ^ nin + 1) 

n=l m=—n ^ ' 



nm I I Hnm \ j 



n + 1 



n[n + 2)(ri — m + l)(n + m + 1) 
(2n+ l)(2n + 3) 



n+l,m) ' 



(20) 



The remainder of the torque is the orbital contribution. 

In optical tweezers, the incident beam that is focussed by the objective 
to produce the trapping beam typically has uniform orbital angular momen- 
tum. In conventional optical trapping, this is a Gaussian beam, with zero 
orbital angular momentum, while if one is intending to rotate trapped parti- 
cles, this could be a Laguerre-Gauss beam carrying non-zero orbital angular 
momentum. In either case, the beam can also carry spin angular momentum 
due to circular polarization. The angular momentum flux of this incident 
beam is determined by the Laguerre-Gauss azimuthal mode and the de- 
gree of circular polarization cXz. The polarization state will, in general, be 
the result of superposition of left- and right-circular modes with az = ±1. 
After focussing, however, the angular momentum is described by the VSWF 
azimutha l mode m. The total a ngular momentum is unchanged by being 
focussed ( Nieminen et al. . 2008bl ). so 



only m = i ±1 VSWF modes can be nonzero. 



for circular polarisation, either the m 
modes will be nonzero, but not both. 



1 modes or the m 



+ 1 



• and both will be nonzero for plane or elliptical polarisation. 

Thus, we will typically have two sets of modes, separated by an angular 
momentum difference of Am = 2. 

While finding the VSWF expansion from the original Laguerre-Gauss 
expansion presents both theoretical and practical difficulties, primarily as a 
result of the VSWF being solutions to the vector Helmholtz equation, while 
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the Laguerre-Gauss modes are solutions to the scalar paraxial wave equation , 
there are a number of ways in which it can be done (INieminen et al.l . l2003bl . 



2007d). 



What then of the effect of the symmetry of the scattering particle? With 
respect to the azimuthal coordinate 0, every object is periodic, even if com- 
pletely devoid of symmetry (in which case, the period is 27r, since the object 
must be rotated by a complete rotation before matching its original ap- 
pearance). Thus, the azimuthal dependence of the fields and the boundary 
conditions (i.e., the effect of the scatterer) present an exact mathematical 
analogy with the case of a planar grating illuminated by a plane wave. For 
the grating, if the component of the wavevector of the incident wave along 
the grating is k^, and the component of the reciprocal lattice vector of the 
grating along the grating is g^, the scattered waves have a discrete angular 
spectrum with wavevectors given by k^j = kx + jqx, where j is the order of 
scattering. This is a fundamental result of Floquet's theorem. In the case of 
the azimuthal behaviour of the fields, the azimuthal mode index m plays the 
part of the x component of the wavevector, and the order of discrete rota- 
tional symmetry p (a scatterer has p-th order discrete rotational symmetry if 
it is indistinguishable p times during a complete rotation) is equivalent to the 
lattice vector. Thus, a discrete azimuthal spectrum, with azimuthal mode 
indices 

rrij = rriQ + jp (21) 

is obtained, where tjiq is the azimuthal mode index of the incident wave. 
Since irij determines the angular momentum of each mode, we see that the 
optical torque is strongly influenced by the symmetry of the scatterer. 

For scatterers that are mirror- symmetric, upward and downward coupling 
must be equal, in the sense that, for example, a mirror-symmetric scatterer 
of 2nd order rotational symmetry (such as a long rod), T-matrix elements 
coupling from m = 1 modes to m = — 1 modes will have the same magnitudes 
as the elements coupling from m = — ltom = l modes. For chiral scatterers, 
these T-matrix elements will, in general, be different. 

Two cases can be simply dealt with. Firstly, if the particle is devoid of 
symmetry (i.e., p = 1), the incident wave is coupled to all possible scattered 
modes, and little can be said about the details of possible effects. Secondly, 
if the particle is axisymmetric {p — >■ oo), we only have rrij = mo, and the 
angular momentum per photon is unchanged. Thus, an axisymmetric particle 
can only experience an optical or electromagnetic torque if the photon flux 
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changes due to absorption (or, in principle, gain). This lack of coupling 
between differ ent angular moraenta is widely exploited in the c a lculation of 
the T- matrix (Nieminen et all 2003a : Mishchenko et al. . 2000l : Waterman . 
197ll : iMishchenkd . Il99ll ). in additions to the benefit of axisymmetry turning 
2D surface integr als into ID l ine in tegrals, or 3D volume problems into 2D 
surface problems (ILoke et al.l . 120071 ). Even in the case of discrete rotational 



symmetry rather than axisymmetry, the restricted coupling between modes 
([2T1) can still be used to reduce computational time by orders of magnitude 
flKahnertl . l2005l : iNieminen et~aD . l2006bl . l2007cf ). 

We can also note some rotation properties of T-matrices. Since the indi- 
vidual VSWF modes have an azimuthal dependence of exp(im0), a rotation 
of the coordinate system of $ about the 2;-axis must induce a phase change 
of m$ in the VSWF modal amplitudes. The combination of rotation of the 
incident and scattered fields is equivalent to a phase change of (mi — m2)$ 
for the T-matrix element coupling modes of azimuthal index mi and m2- If 
the scatterer is rotating, then $ = Qt, and the scattered light experiences 
a rotational Doppler shift of (mi — 1712)^1. This rotational Dopple r shift ac- 



count s for the work done by the field on the spinning particle ([Atkinson 



19351). 



Since the angular dependence of the torque results in alignment to the 
plane of polarization rather than rotation, an object can be rotated at a con- 
trolled rate by rotating the plane of polarization. However, it is not possible 
to control the torque by changing the degree of circular polarization of a 
single beam, since t his alignment will prevent spinning about the beam axis 



( iFriese et al.l . Il998l ). It is possible to overcome this problem by using two 



beams, not mutual l y coh erent, with circular polarization of opposite hand- 
edness ( iFunk et al.l . l2008l ) , but only at the cost of a much more complicated 



trapping apparatus. 

While the use of elongated objects or birefr ingent materials is simple and 



can can allow large torques to be generated (IFriese et al.l . Il998l ). there are 



some drawbacks. Firstly, strongly birefringent materials such as calcite or 
vaterite (both calcium carbonate minerals) can be chemically unsuitable, for 
example perhaps dissolving rapidly in the biological buffer solution being 
used, or not adhering to biomolecules of interest. The shape birefringence of 
elongated objects is typically weak, and the torqu es obtained are typically 



a few percent of h per photon (jBishop et al.l . l2003l ). In addition, elongated 



object tend to align along the beam axis, at which point they are axisymmet- 
ric or close to axisymmetric about the beam, and cannot be rotated further. 
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This last probl em can be overcome by using; flattened o bject rather than 
elongated ones (ICheng et al.l . l2002l : iBayoudh et al.l . l2003l ) , but the problem 
of low torque remains. 



4. J? < 2 versus p > 2: shape birefringence 

If we consider a simple particle such as a mirror-symmetric elongated 
particle, we immediately note that such a particle possesses 2nd order ro- 
tational symmetry, and hence couples modes with mi = m2 ± 2n. We can 
expect mi = 1712 ± 2 to dominate, especially for small particles. 

Firstly, we consider the case of such a particle trapped by a linearly po- 
larized Gaussian beam (m = ±1) mirror symmetric plane polarized beam. 
If the beam is plane polarised along the x-axis, the m = 1 and m = — 1 am- 
plitudes are equal. For simplicity, and without loss of generality, we assume 
that this is the case. If the long axis of the particle, and the symmetry axis of 
the beam coincide, the coupling from m = +1 to m = — 1 must be identical, 
and the scattered modes are of opposing angular momentum are equal. Thus, 
no torque can be exerted on the trapped particle. The same considerations 
apply if the particle is rotated by 90° — in this case, the effect of the rotation 
is a multiplication by —1 of all relevant T-matrix elements, and the scat- 
tered modes of opposing angular momentum are still equal. If we consider a 
single VSWF mode of amplitude a„,i incident on the particle, in the absence 
of absorption, energy must be conserved, and hence Pn,i = an,i(l — ^x) and 
Pn,-i = ix where x is a real number. 

A rotation of 45°, however, results in multiplication by i and — i of the 
upward and downward coupling elements, with no change to the coupling 
from m = 1 to m = 1, or m = —1 to m = —1. As a result, if the initial 
mode amplitudes are real, the real part of, say, the m = 1 mode is increased, 
and the real part of the m = — 1 mode is decreased, while their imaginary 
parts remain the same. Thus, the scattered beam now has non-zero angular 
momentum, and an optical torque is exerted on the particle. 

A left circularly polarized incident beam will simply have the m = 1 
amplitudes reduced, and the m = — 1 amplitudes will become non-zero, again 
giving rise to an optical torque. In this case, the torque is independent of 
the orientation of the particle about the 2;-axis, since the optical angular 
momentum is independent of the phase of the mode amplitudes (which is all 
that will be changed by rotation of the particle). The same conclusion can 
also be simply deduced from the rotational symmetry of the system. 
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In the case of weak scattering, as will be typical of most small elongated 
objects, the coupling will be weak (i.e., x above will be small). For the 
circularly polarized beam, this results in loss of power in m = 1 modes 
proportional to x^, and an equal gain in the m = —1 modes. Therefore, 
the optical torque will be quite small. For a plane polarized beam at 45° 
to the particle, the changes in power in each mode are proportional to x, 
and will generally be much larger. This conclusion is supported both by 
rigorous electromagnetic ca lculations and by ex perimental measurements of 
the torque in the two case (IBishop et al.l . 120031 ). 

Similar considerations apply to the rotation of birefringent particles — 
these possess the same symmetry as elongated particles. The major difference 
is that strongly birefringent particles will results in stronger coupling (i.e., 
larger x above), and the difference noted above between the torque exerted by 
circularly and linearly polarized beams need not apply. Energy conservation 
considerations give the usual maximum angular momentum transfer of 2h 
per photon for a circularly polarized beam, and h per photon for a linearly 
polarized beam. We can note that coupling between the m = 1 and m = — 1 
modes of equal n is not sufficient to reverse the polarization of a beam since 
1 — ix cannot be zero for real x — a waveplate must couple orders of differing 
n in order to be effective. 

If, instead of trapping the elongated particle in a Gaussian beam, we 
now use a beam with an elliptical focal spot, the trapping beam has VSWF 
modes of m = ±1, ±3, ±5, etc, and also modes of higher n, due to the larger 
extent of the focal spot. Again, the overall behaviour will be similar to that 
of an elongated particle trapped in a linearly polarized beam. For weak 
coupling, we can expect higher torques, since the outward coupling to orders 
of greater \m\ now contribute to the torque rather than resulting in equal 
angular momenta in the postive and negative m orders. The exact details, 
however, are likely to be strongly influenced by the exact particle geometry. 
Modes with m = ±1, ±3, ±5, etc will generally exist, due to the overall shape 
of the beam, regardless of the incident polarization, and the interference 
terms that give rise to higher torques in the case of weak coupling will still 
exist — a large increase in the torque should result for circularly polarized 
beams. 

For higher-order rotational symmetry {p > 2), such angle-dependent 
torques are no longer possible. For example, if an object has 4th-order sym- 
metry, clearly we can write the components of its polarizability tensor in the 
plane perpendicular to the beam axis as al, where a is a constant and I is the 
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2x2 identity tensor. Due to the rotational invariance of I, there will be no 
angle-dependent torque. The same argument also applies to any higher-order 
rotational symmetry with p > 2. If the scatterer is also mirror-symmetric, 
there can be no rotation of the plane of polarization of incident light, and 
such a scatterer cannot possess shape birefringence. If the scatterer is chi- 
ral, rotation of the plane of polarization of incident light, and the scatterer 
can have circular birefringence (a.k.a. optical activity). However, circular 
birefringence will not give rise to optical torque. 

This angle-independence of torque can also be seen in terms of VSWFs. 
The two incident circular polarization components are Am = 2 apart in an- 
gular momentum, and a scatterer with p > 2 cannot couple incident modes 
with azimuthal indices m and m + 2 to the same mode. The scattered modes 
resulting from each of the incident polarization components are completely 
independent, and the angle-dependence which results from interefence be- 
tween modes scattered from each of the two polarization components cannot 
occur, since they cannot interfere. 



5. Optically-driven rotors 

One strand of the web of optical rotation has been the practical implemen- 
tation of optical rotation of microstructures, especially specially-fabricated 
microparticles intended for use as optically-driven micromotors, micropumps, 
or other micromachine component s. Work in this field has been recently re- 



viewed by lNieminen et al.l (j2006af ). 

Such optically-driven micromachines typically possess discrete rotational 
symmetry, with p > 2. Therefore, we will consider some qualitative fea- 
tures of the coupling between angular momentum orders. Firstly, we can 
expect scattering to the lowest orders of scattering (i.e., n = 0, ±1) to be 
usually strongest. Thus, most light is likely to be scattered to scattered wave 
azimuthal orders m = mo, mg ± p. 

Secondly, only m such that |m| < nmax are available, with nmax being 
determined by the radius of a sphere required to enclose the micromachine 
element. Therefore, if the order of symmetry p is large, and the microma- 
chine small, only a small number of azimuthal modes will be available. The 
ultimate case of this is the homogeneous isotropic sphere, with p = oo, when 
only the incident m = mo is available for the scattered wave. In principle, 
this can be exploited to maximise torque. As a numerical example, consider 
a structure with p = 8, size such that rimax = 6, illuminated with a beam 
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such that mo = 4. In this instance, the only scattered wave azimuthal orders 
that are available are m = —4, 0, 4, and we might expect such a structure to 
generate more torque, cetera paribus, than a structure that can also scatter to 
m = ±8 (since the scattering to m = +8 will probably be stronger than the 
scattering to m = —8). However, such a structure would have a maximum 
radius of less than a wavelength, which, apart from causing fabrication to be 
difficult, with sub- wavelength resolution being required even at the perime- 
ter, would mean that the device would largely sit within the dark centre of a 
typical tightly focussed optical vortex. Alternative methods of illumination, 
such as a Gaussian or sir nilar beam perpendicu lar to the device symmetry 
axis illuminating one side (IKelemen et all l2006l ). can avoid this, but will not 
necessarily result in greater efficiency. Thus, designing the particle itself to 
more efficiently scatter preferentially into either higher m (i.e., more positive 
or less negative) or lower m (i.e., less positive or more negative) orders is 
desirable. 

Consequently, it is useful to note that, thirdly, the symmetry of scattering 
to positive and negative orders of scattering (i.e., to m < mo and m > 
mo) is dependent of the chirality of the particle shape. If the particle is 
achiral — mirror symmetric about a plane containing the axis of rotational 
symmetry — the coupling from mo = to ±m will be identical, since these 
modes are mirror-images of each other. For a chiral particle, however, the 
coupling to these mirror-image modes will not be identical; even in this case 
of illumination with zero angular momentum modes (i.e., mo = 0), optical 
torque will be generated. If the particle is chiral such that scattering from 
mo = to +m is favoured, then, even for incident modes with mo 7^ 0, 
scattering to higher (i.e., more positive) m will generally be favoured. In this 
way, a particle can be optimised for the generation of torque with a particular 
handedness, leading to greater efficiency when this handedness is the same 
as that of the angular momentum carried by the driving beam. 

Fourthly, we can note another factor that is to some extent a generalisa- 
tion of the second result from symmetry above. Since, for any given particle, 
there are always more modes available with low |m| than with high |m|, 
scattering from a given mo to modes with m = ±5 will usually be such that 
modes with |m| < |mo| will receive more power than those with |m| > |mo|. 
As a result, an arbitrary particle placed in illumination carrying angular 
momentum, but otherwise arbitrary, will usually experience a torque. 

Therefore, we can summarise some general principles for the design of 
optically-driven micromachines: 
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Size matters! The maximum angular momentum available from the 
beam is proportional to the radius of the parti cle, assuming that the 



entire beam can be focussed onto the particle (ICourtial and Padgettl . 



2OOOI ) ■ If the beam is larger, the portion that misses the particle can- 
not contribute to the torque. However, as well as maximising the over- 
lap between the beam and particle, it is also important to maximise 
the angular momentum content of the beam. This can be achieved 
in either of two ways; firstly, by choosing a wavelength such that the 
diffraction-limited Gaussian spot of a circularly polarised beam is the 
same size as the particle (since the angular momentum flux is P/oJ, 
reducing the wavelength maximises the angular momentum), and sec- 
ondly by exploiting orbital angular momentum. The use of orbital 
angular momentum maximises the angular momentum content at a 
particular wavelength, allowing a single set of optical components to be 
used for micromachines of various sizes. We can also note that the vis- 
cous drag on a rotating sphere is proportional to r^. Thus, the highest 
rotation speeds will typically result from smaller micromachines. The 
torque, on the other hand, will increase with increasing size. 

The rotational symmetry of the particle can be chosen so as to optimise 
the torque. The ideal choice depends on the angular momentum of the 
incident beam — for an incident beam of azimuthal order mo, pth-order 
rotational symmetry with p = mo appears to be a good choice, allowing 
coupling of the incident beam to the lowest angular momentum modes 
possible, with m = 0. More generally, mo < p < 2mo should give good 
performance. For smaller p, the difference between the magnitudes of 
the n = ±1 orders is less than 2mo, and thus the difference in coupling 
is likely to be smaller, and for greater p, all scattered orders have a 
magnitude of their angular momentum greater than that of the incident 
beam, which is likely to increase coupling to the non-torque-producing 
zeroth order scattered modes {n = 0). 

Chiral particles can be produced, allowing rotation even by incident 
beams carrying little, or even no, angular momentum. Greater torque 
will also result, compared with a similar particle, when illuminated 
by high angular momentum light of the appropriate handedness. The 
price that is paid is that rotation in one direction is preferred — if one 
desires equal performance in both directions, then an achiral device is 
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necessary. 

• The coupling between incident and scattered VSWFs is essentially a 
vector version of the coupling between an incident paraxial Laguerre- 
Gauss mode and diffracted modes due to a hologram — in both cases, 
the phase variation is of the form exp(im0), and the effect of the sym- 
metry of the particle or hologram on the coupling is the same. This 
allows a simple conceptual model of optically-driven micromachines: 
microhologr ams . 



Figure 1: Simple designs for achiral (left) and chiral (right) optically-driven micromachines, 
based on the design principles outlined here. They can be considered as binary phase 
holograms. The centres are sold for structural integrity. 

The principles outlined above have been applied to the desi gn of optically- 



drive n micromachine elements, and have proved successful (IKnoner et al. 



20071 ). The principles can also be recognised in successful devices that 



appear to be ba s ed on geometric optics principles, or even trial-and-error 



( iNieminen et al.l . l2006al ). Figure [T] shows two very simple p = 4 devices 
based on these principles. The achiral cross rotor is designed to be rotated 
by an incident beam carrying angular momentum, while the chiral rotor 
will rotate in an incident linearly polarised Gaussian beam. Viewing these 
structures as microholographic elements, they can be seen as binary phase 
approximations of interference patterns between paraxial Gaussian and LG04 
beams. If both have planar wavefronts, the achiral rotor results, while if the 
Gaussian beam has curved wavefronts, the chiral rotor results. This suggests 
that the direction of rotation of the chiral rotor can be reversed by reversing 
the curvature of an i nciden t Gaussian beam; this effect has been oberved by 



Galajda and Ormod (120021 ). Simple structures such as those shown in figure 
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[T] have also been fabricated and tested (jUkita and Kanehiral . |2002| ) . The mi- 
crohologram picture of these devices suggests that the optimum thickness is 
that which produces a half-wave phase difference between light that passes 
through the arms of the structure, and light that passes between them. 



6. Optical measurement of optical torque 

One disadvantage resulting from the use of orbital angular momentum 
to drive rotation is t hat while spin angular momentum is relatively easy 



2000 



Nieminen et al.. 2001 



to measure opticallv (iCrichton and Marstonl. 
Bishop et all (20031 . EoO^ iLa Porta and Wane 
cult to optically measure orbital angular momentum. In principle, it should 
be possible to do so by a variet y of methods, such as using holograms 



20041 ). it is rather more diffi- 



as mode filters (IMair et al.l . 1200 ll : iParkin et al.l. 12004). interferometry us- 



ing Dove prisms to introduce a phase shift ( Leach et al.l. 120021) . or naeasure - 
ment of the rotational frequency shift (lAtkinsonl . ll935l : lBasistiy et al.l . l2003l ). 
While such methods have been de r aonst rated to be able to accurately mea- 
sure optical torque (jParkin et al.l . |200J), accurate optical measurement of 



orbital torque within optical tweezers has proved elusive. This appears to 
largely result from sensitivity to alignment, both transverse and angular, and 
aberrations. 

How ever, a method of es timating the orbital torque has been demon- 



strated flParkin et al.l . l2006bl Jah. wherein the rotation rate and spin torque 



are measured for the same object being rotated by left- and right-circularly 
polarised beams, and a linearly polarised beam, all carrying the same or- 
bital angular momentum. If the spin torque is similar in magnitude (the 
handedness will be opposite) for the two circularly polarised cases, and the 
difference in rotation rate as compared with the linearly polarised beam is 
also the same for the two, then we can safely assume that, as a reasonable 
approximation, the orbital torque is the same in all cases. The variation in 
rotation rate with the spin then allows the viscous drag to be measured. The 
orbital torque can then be estimated from the rotation rate in the linearly 
polarised beam. 

It should be recalled that for particles with higher-order rotational sym- 
metry {p > 2), the torques resulting from the incident left and right circularly 
polarized modes are independent of each other. Furthermore, such particles 
do not possess (linear) shape birefringence, and will therefore not alter the po- 
larization of an incident linearly polarized beam. Therefore, the spin torques 
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per photon from each of the left and right circularly polarized components 
of the beam must be equal in magnitude and opposite in direction; if the 
incident beam is linearly polarized, these torques must cancel. 

So, although at first consideration, the above method for the measure- 
ment of the total angular momentum might strike the reader as a hopeful 
approximation, it appears that it is very likely to yield the correct result. 

7. Conclusion 

Prom general principles of symmetry, we have been able to deduce a 
number of features of the generation of optical torque and its dependence on 
the geometries of the particle and the beam. For an elongated or birefringent 
particle, whether rotated by plane polarised light, circularly polarised light, 
or by a beam with an elliptical focal spot, the shared symmetry results in 
similar behaviour in all of these cases. 

Particles of higher-order rotational symmetry, such a microfabricated 
optically-driven rotors ofi^er angle-independent torque, and the possibility of 
higher torques due to exploitation of orbital angular momentum. An impor- 
tant design choice is whether to use a chiral particle, for higher efficiency, or a 
nonchiral particle, for equal performance in either direction of rotation. The 
chiral particle has the advantage of rotating in a Gaussian beam, whereas 
the achiral particle requires a beam carrying orbital angular momentum for 
rotations. In either particle with nth order rotational symmetry is 

best driven by a focussed Laguerre-Gauss mode of helicity of / ~ njl. 

It has been suggested in the past that the essentially unlimited angular 
momentum of optical vortices, with thousands of ft per photon being achiev- 
able, can be used for highly efficient optical rotation. This is not likely to be 
useful for the optical rotation of small objects, since such beams necessarily 
have focal spots hundreds of wavelengths across (if, however, one wishes to 
use an optical rotor that is hundreds of wavelengths across ...). In general, 
the best efficiencies result from the use of focal spots almost the size of the 
particle, with as much power as possible in the higher angular momentum 
modes. 

Pinally, we can note that a randomly-selected irregular particle is almost 
certainly chiral, and should experience an optical torque when illuminated. 
Why, then, is optical torque regarded as unusual, if most randomly-chosen 
particles should experience torques? Although a torque might exist, it might 
also be very small; if smaller than torques associated with the rotational 
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Brownian motion of the particle, the optical torque will not be noticeable. 
Rotation should be more readily observable in environments of low viscosity, 
such as pa r ticles in gas. Under such circumstances, rotation is observed 
(lEhrenhaftl . Il945l ). although convection and thermophoresis may well be the 
dominant effects. Torques on random particle s can also be r aaxim ised by high 
refractive index, increasing the reflectivity (IKhan et al.l . |2006| ). However, 



most naturally occurring light is predominatly unpolarised, and as such can 
be represented as sums of VSWFs with m = ±1. The same is the case for 
linearly polarised light. It is especially in this case that torques can result 
most readily in alignment rather than continuous spinning, which can also 
mask the presence of an optical torque. 
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